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■ Abstract 

■ The Abelian Sandpile Model is a discrete diffusion process defined on graphs (Dhar [10] , Dhar et al. 
which serves as the standard model of self- organized criticality. The transience class of a sandpile 

is defined as the maximum number of particles that can be added without making the system recurrent 
([3]). Using elementary combinatorial arguments and symmetry properties, Babai and Gorodezky (SODA 
^ 2007,[2]) demonstrated a bound of 0(n^°) on the transience class of an n x n grid. This was later improved 

I . by Choure and Vishwanathan (SODA 2012, [7) to O(n^) using techniques based on harmonic functions on 

■ graphs. We improve this bound to 0(n^/logn). We also demonstrate tight bounds on certain resistance 
, ratios over grid networks. The tools used for deriving these bounds may be of independent interest. 

1 Introduction 

J> ' The abelian sandpile model (ASM) is a discrete diffusion process, defined on graphs, pioneered by Dhar [TU] 

to model self-organized criticality in physical systems. It is closely related to the chip firing game investigated 
by Bjorner, Lovasz and Shor [5! and Tardos [24!. This formulation is also known to be closely related to other 

■ interesting phenomena such as stress distribution in earthquakes, size distribution in raindrops, path length 
distributions in loop-erased random walks, for instance. For a nice overview, see the recent comprehensive 
survey article by Dhar The ASM, though easy to define, has a very profound behavior and is far from 

I being completely understood. Research in this area stretches across numerous disciplines such as probability 

theory, algorithmics, theory of computing, combinatorics, non-linear dynamics, fractals, cellular automata, to 
name a few. These connections have been beautifully summarized by Kleber in |16] . Dhar [9J also discusses 
some generalizations of the ASM like the Abelian Distributed Processors (ADP) model which is used to 
model a grid of abstract state machines along with many theoretical and practical applications. 

In the standard sandpile model, "sand particles" are added at the vertices of a (multi)graph. A site 
(vertex) is stable as long as the number of particles at the site remains less than its degree. Adding more 
particles would render the site unstable and is accompanied by the unstable site's passing a particle along 
each edge to its neighboring sites. This relaxation process is referred to as toppling. One of the sites known as 
the sink cannot topple. To ensure that every relaxation process eventually stabilizes, one needs the condition 
that the sink is reachable from every other site. As the system evolves, the sandpile goes through a a sequence 
of configurations. Those which can be revisited in any toppling sequence are called recurrent^ the remaining 
ones are termed transient. Typically, one starts with the empty configuration and as particles are added, one 
moves through transient configurations till a recurrent configuration is reached. Thereafter the configurations 
stay recurrent. The steady state behavior of a sandpile is characterized by its set of recurrent states. It has 
been observed by physicists that for most natural phenomena, the time taken to reach a recurrent state is 
small. Hence any acceptable model must refiect this tendency to reach steady state rapidly and it becomes 
important to study the time taken to reach recurrence in these models. 

The essential parameter in our discussion is the number of particles which ensure recurrence. If particles 
are added randomly then a simple coupon collector type argument demonstrates polynomial bounds on the 
expected time to recurrence (as already mentioned in [2J). The other scenario is to add particles adversarially 
so as to avoid a recurrent state for as long as possible. This problem was highlighted by Babai and Toumpakari 
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[3] where they define the requisite number of particles as the transience class of the sandpile. This later 
motivated the insightful work by Babai and Gorodezky [2\ on grid based sandpiles which are the most 
studied objects as compared to any other graph class because of their outstanding significance in statistical 
physics. In their path-breaking paper, Babai and Gorodezky [2\ show that for the standard n x n square 
grid based sandpile, the maximum number of particles one can add before hitting a recurrent state is 0{n^^). 
This is a remarkable result in view of the fact that some closely related sandpiles (for example line graph 
based) have transient state paths of length exponential in graph size. They use intricate combinatorial 
arguments based on particle conservation and the symmetry group of grid graphs to demonstrate the above 
mentioned bounds. Choure and Vishwanathan [7^ obtained general bounds on the transience classes in terms 
of harmonic functions on graphs. They improve the upper bounds to O(n^) and obtain the first non-trivial 
lower bounds of il{n^). They also demonstrate a very general bound on the transience class in terms of 
sandpile size and the number of connections to the sink and explicit algebraic expressions which bound the 
transience class values of planar sandpile. However, simulations suggest a bound close to O(n^) for the grid 
sandpile. Also, it made apparent that every improvement in the current known bounds would culminate in 
better understanding of the model. 

Our contribution: We build upon the results in [7 and improve the bounds by a logarithmic factor to 
0(n^/logn). This is achieved by using fundamental results from classical (electric) network theory, along 
with a bit of geometric graph theory. For the grid circuit corresponding to the grid sandpile, we obtain 
tight bounds on the effective resistance across several pairs of vertices. This is used to bound the resistance 
eccentricity of the grid network and improves upon the previous known bound by a logarithmic factor. Finally, 
using the basic results from [7J, this improvement implies an improvement in the transience class bound. We 
believe that the techniques used by us will be of independent interest. This marks a departure from most of 
the prevailing work on bounding the effective resistances which relies on algebraic tools or results based on 
random walks. 

Related Work The classical theory of electric networks along with the well understood connections 
with random walks ([20], [12], [17]) has some very powerful and intuitive results. These results have recently 
found applications in almost every important area of theoretical computer science |6], |15j . [23], [T]. [^ use 
the theory of harmonic functions to analyze sandpile behavior (in the context of diffusion) in analogy with 
the theory of random walks on graphs. The results we report in this paper are an extension of their work 
and achieves better bounds using some powerful theorems from electric network theory. Advances with a 
complexity theoretic flavor include proof of the one-dimensional sandpiles prediction problem in LOGDCFL 
by Peter Bro Milterson [TS]. Also. Schulz [2T] mentions a related NP-complete problem. The group structure 
of the space of recurrent configurations, first introduced by Dhar. Ruelle, Sen and Verma in [llj, is also a 
fertile area of analysis. Cori and Rossin [5] show that sandpile groups of dual planar graphs are isomorphic. 
Toumpakari [26] discusses some interesting properties of sandpile groups of regular trees where questions 
related to group rank are studied and the paper is concluded with an interesting conjecture on the rank of 
all Sylow subgroups of the sandpile group. Specific families of graphs like square cycles C^, x C„, 3 x n 
twisted bracelets, etc have been analyzed in [14], [22], |13| . 

2 Preliminaries 

2.1 Basics of Harmonic Functions and Potential Theory 

For a very nice introduction to harmonic functions on graphs, we refer the reader to the beautifully written 
paper by Benjamini and Lovasz |3] and to Teles j25J for a thorough view. We start with some important 
definitions and fundamental properties. Given a connected graph G and a function tt : V{G) M, we say 
that TT is harmonic over Vh if. 



The remaining vertices (lying in V — Vh) are called the "poles" of tt. The set Vh is also called the interior 
of TT with vertices adjacent to the set of poles referred to as the boundary. We see that the value of tt at any 
vertex in Vh is the average of its value in the immediate neighborhood. In case of multi graphs, we take the 
appropriate weighted means, where the weights are the number of common edges. This leads us to the first 
basic property. 
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Property 2.1. Any non-constant harmonic function can assume its extreme values only at the set of poles. 

It follows that every non-constant harmonic function has at least two poles, its maxima and minima. 
Such functions are completely determined by their values on these vertices. Formally speaking, 

Property 2.2. Uniqueness: If two functions harmonic on Vh agree on the boundary, they agree everywhere 
in the interior. 

More generally, we have the following property. 

Property 2.3. Given a set of poles, a harmonic function is uniquely determined modulo scaling and trans- 
lation by a constant. 

Properties 12.21 and 12.31 important as they allow one considerable freedom in constructing harmonic com- 
pletions of functions defined over the boundary set. This problem is the discrete analogue of the classical 
boundary value problems in complex analysis. We will describe two important examples in which these 
function arise naturally. 

Random Walks on Graphs: Consider a graph G and two special vertices s and t. The potential associated 
with V, with s and t as poles, sT^tiv) is defined as the probability of reaching t before s starting from v. 
One can check that the function tt so defined is indeed harmonic on the set V — {s,i}, with the maximum 
value of 1 at the node t and the minimum value at s. The generalization to the multi-pole situation is also 
straightforward. 

Electric Networks: Consider a resistive electric network (i.e. a circuit made up entirely of resistors). Let 
s7rt(f) be the potential that appears at node v when unit potential is applied across t and s. Using the 
equation of charge conservation (KirchofF's node law), one can show that these potentials are harmonic on 
all nodes except s and t. 

The main implication here is that one can intuitively think of the electric network theory as an analysis 
of random walks of electrons on the underlying graphs. Consequently, results from network theory can be 
used to prove interesting facts in other related areas. As an example, consider the problem of constructing 
the harmonic completion of a function with given boundary values. All one needs to do is to take the 
corresponding circuit and apply potentials equal to the boundary values on the boundary points. The 
potentials that will appear on other nodes can be computed using basic linear algebra (the only non-trivial 
step involves inverting the combinatorial Laplacian of G) thus allowing construction of harmonic completions 
efficiently. We present below two very basic results on network analysis which will be needed in the following 
sections. 

The effective resistance between a pair of nodes u and u, Reff{u,v) is defined as the potential difference 
which develops between u and w if a unit current source is applied across u and v. 

Lemma 2.1. Potential Reciprocity Lemma : If taking s and t as poles with 7r(s) = and 7r(i) = 1 induces 
a potential of gTit{v) at node v and interchanging the roles of v and t induces s'^v{t) o,t t then, 

Reffis,t)sTrt{v) = Reff{s,v)sny{t) (2) 

Proof : Consider the given network G with the special node s G V{G). We refer to the corresponding 
modified network G(e) obtained from G by adding an edge with resistance 1/e between every node and s. In 
particular, G(0) = G. Furthermore, we refer to an edge between s and u by sli. We will be using the current 
source version of the reciprocity theorem. If applying a unit current source across st results in a potential 
of V across sv, then applying unit current source across sv results in a potential of v units across st. The 
value of this potential v can be expressed, using Ohm's law, as the ratio of current through the edge and the 
resistance of the e— edge between the particular node and sink. Since both potentials are equal in magnitude, 
we can say that on G(e), 

Reff{s,t)sTTt{v) = Reff{s,v)s'^v{t) 

This follows from observing that applying a unit current source across sv is equivalent to applying a 
voltage source of R(,ff{s,v) across sv. Because of linearity, it follows that a potential of s''^^{t)Reff{s,v). 
appears at node t. Similarly so for the other configuration. 

This equation holds for arbitrarily small values of e. Consequently it holds for graph G(0). ■ 
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In particular, when the effective resistances across s and t are the same as s and v, we have s7rt(w) =s T^vit)- 
In the following discussion, we will omit the left subscript (s) from ^TTt whenever it is clear from context. We 
say that a walk P is an instance of gT^t if it starts at some vertex w, avoids s and ends at t. For a proof of 
the following lemma, the reader is referred to [7]. 

Lemma 2.2. A triangle inequality for potentials 

TTi[j).TTj{k) < TTi{k) (3) 

Remark: The utility of this inequality becomes clear when interpreted in the context of electric networks. 
Consider a network such that the node with ground potential is fixed and we are allowed to apply power at 
any other node and observe the resulting potentials. The inequality implies that if applying a potential Vi 
at i produces unit potential at node j and applying V2 at node j produces unit potential at node k. then 
applying V1.V2 units at i produces at least unit potential at node k. 

2.2 Introduction to the Abelian Sandpile Model 

Our notation and terminology follows Babai and Gorodezky [2]. 

Definition 1. A graph G is an ordered pair (y(G'), E{G)) where V{G) is called the set of vertices and E{G) 
is a set of 2— subsets of V, possibly with repeated elements, the set of edges. 

This is referred to as a multi-graph in literature but we will use graph for brevity. The degree of a vertex 
V V is defined as the number of edges in E which contain v. Two vertices v and u are called adjacent (or 
neighboring) if (w, v) G E. A path between two vertices u and v is an ordered sequence of edges ei, 62, . . . , Cfc 
such that u d ei, V € Ck and for all values of i, e; H e^-i-i ^ <j). The graph G is connected if there exists a path 
between any pair of vertices. 

To model an Abelian Sandpile Model, we take a connected graph G with a special vertex called the 
smfc, denoted s G V . Non-sink vertices in G are called ordinary vertices and this subset will be denoted by 
Vo = V- {s}. 

Definition 2. The configuration of a sandpile G is a map c : Vo ~> N, which will be represented as a vector. 
The weight of c is |c| = Y,vev„ "^(^)- 

The configuration c records the number of sand particles contained in each of the ordinary sites. The 
empty configuration is the zero vector. The capacity of a site is the maximum number of particles that it can 
hold and is one less then the degree of the node. 

Definition 3. An ordinary node v is said to be unstable in a configuration c if c(v) > degree(v). The 
configuration c is said to be unstable if any site under it is unstable, else it is referred to as stable. 

When a site is unstable it is said to topple, that is it passes on some of its particles to its neighbors. 
When a site v topples once, it loses degree(v) particles and each neighbor of v acquires a particle for every 
edge common with v. The sink node never topples. Starting with the empty configuration, we keep adding 
particles one by one on sites of our choice and topple them when necessary. 

The ASM evolves in time through two modes, particle addition at sites and relaxation of unstable sites via 
topplings. A toppling sequence is an ordered set of configurations where every configuration can be obtained 
from the previous one by toppling some unstable site in it. Note that the event of many sites becoming 
unstable simultaneously poses no complication since the order in which they are subsequently relaxed does 
not affect the final stable configuration that is obtained at the end of toppling sequence. Elementary proofs 
of such confluence properties can be found in the pioneering paper on ASMs by Dhar |10| . See also Babai 
and Toumpakari [3J. 

A configuration is called recurrent if it is reachable from any configuration. As already mentioned, we 
say that a configuration Ci is reachable from a configuration Cj if by adding some particles to Cj (possibly 
at multiple sites) and subsequently relaxing it, we can obtain c^. A configuration is transient if it is not 
recurrent. The set of recurrent configurations is therefore, closed under reachability. 

Property 2.4. // 3c' such that there is a toppling seqquence from c' to c in which every site has toppled at 
least once, then c is recurrent. 
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The proof follows from the fact that the existence of such a toppling sequence precludes the existence 
of forbidden sub-configurations and hence makes the configuration recurrent. For a complete discussion on 
forbidden sub-configurations and recurrence of configurations the reader is referred to [10^ and [9J . 

We analyze the process of adding one grain at a time to the sandpile and study its evolution. As in 
the standard theory of Markov chains, recurrence characterizes the long term (steady state) behavior of 
sandpile. Our investigation is concerned with the maximum number of particles that can be added while 
staying transient. Following Babai and Gorodezky [5], for a sandpile S we define, 

Definition 4. The transience class of S denoted by tcl(S'), is defined as the maximum number of particles 
that can be added to S before reaching a recurrent configuration. 

In view of propertv l2.41 we can bound the transience class from above by the maximum number of particles 
that can be added before all the nodes have toppled at least once. 

For the sandpile S, consider the corresponding circuit obtained by taking the resistive network formed 
by the same underlying graph with a unit resistance along each edge. Also, the potential of the sink node 
is fixed at zero. The potential at any node v when unit potential is applied at node w is denoted by iTw{v). 
Choure and Vishwanathan [7] prove the following bounds using results from LP duality, network theory and 
grid symmetry. 

Theorem 2.1. (J^) tcl{S) is 0(max„,^,{7r^(^;)-i T,vWv) - l).7r^{v)}). 

This in turn implies the following weaker version. 

Corollary 2.1. tcl{S) is 0{ma,Xy^^{Tr.u]{vy^}.ma.Xy^^{J2y{d{v) - l).Tr^{v)}). 

This allows us to obtain separate bounds on the contributing factors and use them to obtain (possibly 
weaker) bounds on the transience class. [7] also show that, while finding worst case bounds, when it is enough 
to consider only the nodes directly connected to the sink. From the following sections, we will now continue 
this discussion in the context of grid sandpile and assume that both v and u are boundary nodes on the grid. 

3 The case of Grid Sandpile 

As noted in the introduction, the sandpile associated with the n x n grid is of particular importance. We 
define it formally below. 

Definition 5. Consider the n x n grid graph. Attach an extra sink node to the boundary such that there is 
a single edge to each non-corner boundary node and double edges to the corner nodes. We denote both the 
sandpile and the corresponding circuit by GRID„. 

Notation:Nodes on the grid are labeled with variables u,v,w, . . .. The sink node is labeled s and the node 
at the center of grid (if it exists) is called c. For an illustration of the sandpile and the corresponding resistive 
network, see figure ([1]). 



All edges have 
one ohm resistance 



'light edges connect 
to the node s 



Figure 1: The graph of GRID„ 

Note: In the following discussion, we will be assuming that n is odd, ensuring the existence of a center site. 
It is however easy to extend the discussion to the case of even values of n. 
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Babai and Gorodezky [2] have shown that id(GRID„) = 0{n'^^ °°^''). This was later improved by Choure 
and Vishwanathan [7] to 0{n'^) using the theorem (|2.ip . In this paper we will improve this bound to 
0{ii7 /log{n)). As mentioned earlier, the bound is product of two components : riu(GRID„) = '^^{d{v) — 
l).'Kw{v), the potential profile, and the maximum value of t:t,{w)^^ , over all (boundary) values of v and w. 
The following lemma from [7] bounds the potential profile. 

Lemma 3.1. (I7j) Consider any vertex v on the boundary of grid in GRIDn- The potential profile Tw{GRIDn) 
induced due to 0(1) current source at w is 0{n). 

The maximum value of 7r„(z/;)~^ is that value of potential, which when applied anywhere on the grid, 
induces at least unit potential on every other node. The following lemma in [7J bounds this value. 

Lemma 3.2. (I7l)ln GRIDn, applying 0{n^) potential at any site induces at least 0(1) potential everywhere. 

The last two lemmas, along with theorem (|2.ip yields the O(n^) bound of [7J. We will improve this by 
improving the bounds on max7rt,(?x;)~^. 

3.1 Tightening the max7r„(iy)^ bound 

We will now obtain a lower bound on the minimum value of 7rt,(z«) for any pair v and w occurring on boundary. 
Let the center of the grid be denoted by c. Then, using Lemma 12.21 (triangle inequality of potentials), we 
obtain nw{c)T:c{v) < t:v{w). Using Lemma |2.H 



T^wic) = n^iw).—-^ 

Reff(s,w) 

Clearly, /Jmin^, ndv)'^ < min„^^ Tr^{w), where /? is the minimum, up to constant factors, value of 
over all possibilities of boundary nodes w. 

Lemma 3.3. (171) If applying the potential p(n) on a corner induces unit potential at the center, applying 
K.p{nY I j3 at any node induces unit potential at every non-sink node, where K is a constant. 

The following lemma from [7J bounds the value of p(ri). 

Lemma 3.4. In GRIDn, applying 0{n^) potential at (1,1) induces at least 0(1) potential at the center. 

[7] demonstrate a constant lower bound on the value of /3, which results in an 0{n^) bound on max7r„(?x;)~^. 
To achieve tighter bounds on the transience class of grid sandpile, we will show that /? is in fact 17 (log n). 
We start with mentioning the following fundamental principal in network theory. 

Theorem 3.1. Rayleigh's Monotonicity Principal: In a resistive network G, consider any pair of 
nodes u and v. The effective resistance between them increases monotonically as a function of any particular 
resistance. 

In particular, if a given pair of nodes are fused into one, the effective resistance across any pair can not 
increase. This is because, the fusion operation can be seen as reducing the edge resistance between any two 
nodes to zero. Similarly, cutting open any edge resistance in equivalent to increasing that resistance value 
to infinite. This operation cannot decrease any effective resistance, the following corollary summarizes this 
simple observation. 

Corollary 3.1. In a resistive network G, consider any pair of nodes u and v. The effective resistance between 
them cannot increase is any two nodes in G are fused and cannot decrease if any particular edge is deleted. 

To bound j3. we separately bound the value of R^f f{s,w) from above and Reff{s,c) from below. An 
upper bound on R^f f{s,w) is easy to obtain. 

Lemma 3.5. R^f f{s,w) is bounded above by unit resistance. 

Proof :The constant upper bound follows from the fact that it is a parallel combination of a unit resistance 
with the remaining network. The net resistance of parallel combination of ri and 2 is at most min{ri,r2}. 
Therefore, Reff{s,w) < 1. □ 
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Lemma 3.6. Reff{s,c) is n{\ogn). 

Proof; This proof will involve repeated application of corollary 13.11 till we obtain a circuit in which 
effective resistance is easy to bound. We start with the network shown in figure [TJ In this network, there are 
two axis of symmetry parallel to some edge, as shown in figure [51 For any node v not on any of these axis, 
there exist three more which can be obtained via (possibly multiple) reflections in these axis. The first step 
is to perform a merging of these four nodes, for each node in the original grid. For each of the nodes lying on 
the axis, there is just one more counterpart, and it lies on the same axis. Fuse this pair, for all nodes lying 
on the axis. 



Figure 2: GRID„, identifying all four symmetric points, leads to folding along the axis of symmetry shown 

This gives us a circuit which looks like the one figure [3l However, the edges on the left and top have half 
units of resistance each while all others have a quarter unit. This is because the edges on the top and left 
are formed by a parallel combination of just two edges while all others come from a combination of four. To 
remedy this, we apply the corollarv l3. II and reduce these edge resistances to quarter units without increasing 
any effective resistance. This gives use the circuit in figure [3l 



all edges have 
resistance of .25 units 



Figure 3: Circuit obtained after identification, sink is attached to only two of the sides 



The next operation also involves fusing but of a different type. Refer to figure We will now be fusing all 
vertices which are at the same shortest path distance from c. Hence, all vertices with distance d are merged to 
form new vertex, Vd- One can also view this as reducing the resistance values between some special diagonally 
opposite vertices, from infinity to zero. These resistances are shown in red in figure 2] 

This fusion gives us the circuit shown in middle. For aid of illustration, let us number the resistances 
starting from c. So, ri is the resistance between c and wi, r2 between vi and V2- In general, the resistance 
Ti connects Vi-i and Vi. Evidently, each resistance ri{i < n/2) is formed from a parallel combination of 2i 
resistors each of value of 1/4 units. This gives n = l/2i.K, where K = 1/4. The figure shows the resistance 
values in multiples of K. For values i > n/2, it is straightforward to work out similar bounds. 

The final step involves fusing all the vertices directly adjacent to sink. This yields the last circuit in the 
sequence shown in figure It consists of the node c connected to sink node through a series combination of 
ri with i ranging from 1 to n/2. 
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C 1/2 1/4 1/6 ... 1/n S 

Figure 4: Reducing the grid like circuit to a simple line circuit 
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Following Reiligh's principle, the value R^ff forms a lower bound on the value of Reff{s,c). The lemma 
follows. □ 

Using the bounds shown in Lemmas 13.51 and 13. 6[ we obtain the following bounds on l3. 
Lemma 3.7. The value of (3 defined above is lower bounded by Q([og{n)). 

Remark: It is evident that the reductions performed above to derive the lower bounds are some what liberal 
and one might be inclined to believe that by a more careful analysis, these bounds can be improved. However 
this is not so. Using the monotonicity principal, one an cut the appropriate edges so that the remaining 
configuration has an n/2 x n/2 grid with the center and sink attached to diagonally opposite ends. This 
configuration ensures that the effective resistance is bounded from above by 0(log(rt)). The logarithmic 
upper bounds are folklore in random walk theory, see [19] for details of obtaining these bounds. This shows 
that /3 is in fact 8(log(n)). 

3.2 Transience Class of GRID„: a new bound 

Using Lemmas l3.3| 13.41 and 13. 7| we obtain the following result which bounds the value of maxi,_u,(7ri,(w))^^ 
from above. 

Lemma 3.8. In GRIDn, applying 0(n^/ log n) potential at any site induces at least 0(1) potential everywhere. 



Using Lemma [3. 8113. H and Theorem 12. II we have the following bounds on teZ(GRID„). 
Theorem 3.2. td{GRIDn) = 0(n^/logri). 

4 Conclusion 

The main open question is that of tightening the bounds on tc/(GRID„). The general prevailing perception 
is that one can expect substantial improvements only in the estimation of max„^„ tTj^{v)~^. We believe the 
approach using the triangle inequality of potentials which uses the grid center as potential pivot is close to 
its limits. The improvement can come through probably constructing a better harmonic distribution on grid. 

Another important direction is showing better lower bounds. [7] demonstrate that the transience class of 
grid sandpile is il{n^). Following empirical evidence, the true value is going to be around n{n'^~'^) where e 
is a small positive constant. The main reason for the gap between the observed and known bounds is that 
no non-trivial lower bounds on the potential profile are known. It seems very likely that the potential profile 
would obey bounds of the type n(n^~'^). These would be enough to imply the required lower bounds on the 
td(GRID„). We refer the reader to [2] [7] for a fuller discussion on various other open problems. 
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